METHOD OF DETERMINING THE VOLUME SCATTERING FUNCTION OF 
OCEAN WATERS IN THE BACKWARD DIRECTION USING A SATELLITE 

OCEAN COLOR SENSOR 

Field of the Invention: 

The present invention provides a method and a system for determining Volume 
Scattering Function in rearward direction of ocean waters using a Satellite Ocean Color 
Sensor The present invention also provides a method and a sytem for determining 
spectral components of the volume scattering function for different types of ocean water 
and thereafter recognize these waters by their optical scattering properties. 

Background and Prior Art References to the Invention: 

The remote sensing of ocean color by use of a radiometric sensor on a satellite platform is 
a widely used technique to determine the concentrations of major seawater constituents in 
15 the ocean. These constituents are phytoplankton and their detritus, inorganic suspended 
material and dissolved organic matter. A remote spectral radiometer on a satellite 
platform is tuned to receive light at several wavelengths in the visible region and near- 
infrared portions of the electromagnetic spectrum. They are several pathways by which 
photons from the sun can enter the 'field-of-view ' (FOV) of the satellite sensor after 
20 interacting with the ocean surface and the atmosphere. 

The pathways are illustrated in Fig.l. For example path (a) illustrates direct specular 
reflection of the sun's rays from the ocean surface; path (b) corresponds to direct 
scattering from the atmosphere alone, while path (c) represents rays that are scattered 

25 both by the ocean surface and by the atmosphere. The route followed in dashed path (d) is 
that taken by upwelling light backscattered by different oceanic layers that ultimately 
reach the ocean surface [2] and are further attenuated by the intervening atmosphere 
before arriving at the satellite detector [3]. The radiance component (d) carries useful 
information on the optical properties that are related to seawater constituents that color 

30 the ocean. 

Most satellite sensors have the means of orienting the detector so as to avoid or minimize 
specular reflection of path (a) The atmospheric signal arising from path (b) constitutes 
nearly 80% of the signal scattered into the satellite detector, and carries no useful 
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information about the ocean. There are techniques which have been developed by Gordon 
&Wang (1994), Gordon (1997), and Gordon & Morel (1983) which deal with the 
removal of the atmospheric component. This will not be of concern to us here. 

5 Hitherto known devices to measure general-angle scattering meter was by Tyler and 
Richardson (J.E. Tyler , W.H. Richardson, " Nephelometer for the measurement of 
volume scattering function in-situ", J. Opt.Soc. Am 48, 354-357 (1958)). The first fixed 
angle optical backscatter devices were reported by Moore et al (C.A. Moore, R.Honey, D. 
hancock, S.Damron, and R.Hilbers, Development and use of computerized optical sea- 

10 truth instrumentation for LIDEX-82, SRI International Project 3878, Final Rep. 1984) , 
and substantial improvements on these devices undertaken by R. Maffione and D.R. 
Dana (R.A. Maffione , D.R. Dana , " Instruments and methods for measuring the 
backward-scattering coefficient of ocean waters Applied Optics, Vol.36, No. 24 pp 
6057-6067, 1997) with the purpose of finding the best estimate of the back scattering 

15 coefficient (b b ) via a measurement of Volume Scattering Function (VSF) at a fixed 
backscatter angle of 142°. The optical geometry of the backscatter meter developed by 
Maffione and Dana is illustrated in Fig.2. The meter consists of built in collimated 
source of light [a light emitting diode - LED] having a finite divergence and propagated 
at a fixed refracted angle through a glass window of the instrument casing into 

20 surrounding seawater. The sample volume at which backscattering occurs (shown as a 
shaded area in fig.2) is the geometric intersection of the field-of-view [FOV] of the 
source and detector. The detector is placed adjacent to the source. The backscattered flux 
that falls within the solid angle cone subtended by the sample volume at the detector will 
be the signal measured by the detector. The receiver detector is frequency modulated so 

25 that it phase locks to only the backscattered flux from the LED source. This electro- 
optical technique rejects ambient light changes and permits its use in bright sunlight. The 
backscatter coefficient bb in this meter is computed on the basis of the Integral Mean 
value Theorem such that there is a fixed scattering angle \\j * of the VSF function where 
b b = 27t.p (i[/ *). The detected signal is proportional to P (\\f *) multiplied by a weighting 

30 function W (z;c) which varies with the distance z of the sample volume from the 
backscatter window (see fig.2), and the attenuation coefficient of seawater. Most 
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backscatter devices are deployed from a ship or a boat. There are several known 
drawbacks to these optical backscatter devices which are identifiable 

• The small volume from which the backscattered flux originates may not capture a 
representative distribution of suspended particles in the ocean. 

5 • The detected signal has to calibrated by a separate experiment using a Lambertian 
target in clean water tank. 

• The beam attenuation coefficient [c] of the unknown water needs to be known so as 
to apply a correction to the measurement of backscattered flux. Further if [c] changes 
by a large factor, as happens in coastal waters, then the centroid angle (\\f *) at which 

10 the volume scattering function is measured will also change by several degrees. 

The present invention recognizes that the sun, the ocean and a satellite ocean color 
sensor constitute the elements of a giant backscatter sensor system in the sky. The source 
of radiation is the sun ( the analogue to this is the light emitting diode as the exciting 
15 source of the backscatter meter), the ocean surface is the sample volume which emits 
backscattered flux from the top layer and all layers below the ocean surface in response 
to the stimulation provided by the sun's radiation, and the detector in our system is the 
satellite ocean color sensor which receives the backscattered flux from the ocean ( the 
analogue to this is the single photodiode detector used in the backscatter meter ). 

20 

The present invention provides the means to interpret the satellite signal in terms of the 
Volume Scattering Function [VSF] in the backward direction at fixed angles by using the 
existing geometry of a known sun-ocean -satellite detector system. Our best example is 
the Ocean Colour Monitor (OCM) on the Indian satellite platform IRS-P4, although the 

25 same ideas contained here could be applied to other Ocean Color Sensors namely 
SeaWiFs on the SeaStar satellite. The OCM operates as a push-broom camera which 
uses as detector, a charge coupled linear array device (CCD) as a light detector (see 
Fig.3). The various pathways of photons (as depicted in Fig.l) result from reflected and 
backscattered events in the atmosphere and from within the ocean's interior. These 

30 pathways are contained within the cone defined by a field of view of + 43 °, and a total 
swath width of - 1420 kms. at sea level as seen by the satellite optics ( Fig.3 again ). 
However, within this large cone, each element of the CCD array views an ocean pixel 
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with infinitesimally narrow cones having fixed look angles and a IGFOV (Instantaneous 
Ground Field of View ) of 360m x 236 m. There are 3730 active elements on the CCD 
array viewing an equal number of contiguous ocean pixels spanning the swath width of 
the satellite track. Each scan line in the push broom operation is contiguous with each 
5 other. Table 1 summarizes the major features of the OCM sensor as it exists on IRS-P4. 
The eight spectral channels on OCM have separate collecting optics and separate CCD 
line array detectors each having their own drive and signal processing electronics. 



Table 1 : Features of OCM Camera ( see reference 5 ) 



Feature 


Description 


Orbit 


Polar sun synchronous 


Satellite Altitude 


720 kms 


Swath 


>1420 kms 


Repeativity 


2 days 


IGFOV (Ocean pixel) 


360m x 236 m 


Spectral Channels (nm) 
Bands 1 -8 


414.2, 441.4, 485.7, 510.6,556.4, 669, 
768.6, 865.1 


Active CCD pixels 


3730 elements 
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Objects of the Present Invention: 

The main object of the present invention is to provide a method for determining the 
volume scattering function of ocean waters in the backward direction using a satellite 
ocean color sensor thus enabling large scale coverage of the scattering. 

is 

Another object of the present invention is to demonstrate that geometric anlysis of OCM 
sensor optics results in an operational algorithm that uses existing satellite products to 
provide a new way to visualize optical scattering of the ocean surface. 

20 Still another object of the present invention is to demonstrate by example the variation of 
the VSF function at 490nm using real satellite data 
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Yet another object of the present invention is to show that the form of the operational 
VSF function can be extended to other visible wavebands of the sateliite ocean colour 
sensor provided valid data on the other satellite products are available. 

5 In a further object of the present invention, the method provided here could be applied to 
other Ocean Color sensors as has been done in the data presented in fig.6 using the 
SeaWiFs sensor. 

Brief Description of the Accompanying Drawings: 

10 In the drawings accompanying this specification, 

Figure 1 illustrates the different paths taken by the sun's rays in reaching a satellite 
detector. This has been describes earlier in the opening paragraph. 

Figure 2 represents the optical geometry of a commercial backscatter instrument . The 
15 geometric intersection of the field of view (FOV) of the source ( light emitting diode ) 
and the detector shown shaded is the sampling volume from which backscatter radiation 
emanates. The source and detector are both mounted on an endcap of the instrument. 

Figure 3 represents the satellite sensor system [3] looking down at the ocean surface [2]. 

20 The Instantaneous Field of View ( IFOV) of an ocean pixel [4] is 360 meters across track 
and 250 meters along track. The swath W of the satellite is 1420 kms, and the altitude H 
of the satellite is 720 kms. This produces a cone of half angle (0 ) equal to 43 degrees 
subtended by the swath W at the CCD array detector of the satellite sensor. Individual 
pixels on the CCD array view water pixels of the above size with infinitesimally small 

25 FOV's having near zero divergence angles. 

Figure 4 represents the view seen by a center pixel of the satellite's CCD detector of the 
looking straight down at an ocean pixel of size 360m x 250m. The angle of the sun's rays 
is represented by (0 S ) and the angle of the upwelled flux (*) to the zenith direction that 
30 reaches the satellite pixel is shown as (0 P ). The backscattering angle (\(/) is the 
deviation produced in the original solar rays. Note that the upwelled flux is an integrated 
flux composed of all backscatterd events from different layers in the upper ocean under 
solar excitation. 
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Figure 5 shows a high resolution variation of the VSF function at 490nm as a function of 
the backscattering angle along a scan line of an OCM image of 16 March 2002 where 
Trichodesmium blooms were viewed in abundance. 

5 Figure 6 shows plots ( upper) of the backscatter coefficient at 490nm measured at 
different stations on a research cruise in November 2001. The corresponding VSF 
function at 490nm ( lower plot) was computed using SeaWifs data obtained from 
matched overhead passes at each of the stations. 

10 Summary of the Invention: 

Accordingly, the present invention provides a method of determining the Volume 
Scattering Function of Ocean Waters in the backward direction using a satellite Ocean 
Color Sensor. The said method exploits the geometry of the sun-ocean-satellite detector 
to function as a backscatter system comprising of the sun [1] as a source of radiation, the 
15 ocean as the sampling volume [2], and the satellite ocean color sensor [3] as the detector 
of the backscattered flux from the ocean emanating from separate ocean pixels arranged 
on a scan line across the track of the satellite payload and arriving at corresponding 
pixels on the CCD array detector of the satellite sensor [3]. 

20 Detailed Description of the Present Invention: 

Accordingly, the present invention provides a method for determining the Volume 
Scattering Function (VSF) of ocean waters in backward direction using a satellite ocean 
color sensor, said method comprising the steps of: 

(a) obtaining satellite view angle (9 sat ), solar angle (0 s ), remote sensed reflectance 
25 [Rrs ( 0 P , X)] and diffuse attenuation coefficient (K^) at each ocean pixel on a 

scaline using a satellite ocean color sensor; 

(b) calculating corrected look angle of a pixel of satellite ocean color sensor on an 
ocean pixel using: 

0 P = Qsat + sin _1 (0.113 tan 0 sat ); and 
30 (c) calculating the Volume Scattering Function of ocean waters in the backward 
direction at a given wavelength [p(v|/, A,)] using: 

P(i|/, >,) = Rrs(0 p »Kd [cos 0 S cos0 p ] 
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In an embodiment of the present invention, the satellite ocean color sensors include CCD 
array detectors, ocean color monitors and SeaWifs Oocean color sensors. 

In another embodiment of the present invention, the detectors view a scan line of ocean 
5 pixels arranged along a swath of the ocean surface. 

In yet another embodiment of the present invention, the swath "W" of the satellite is 
about 1420 Kms. 

10 In still another embodiment of the present invention, about 3730 ocean pixels are covered 
in a satellite scan line. 

In one more embodiment of the present invention, wherein Instantaneous Field of View 
(IFOV) of an ocean pixel about 360 meters across track and about 250 meters along 
15 track. 

In one another embodiment of the present invention, the Volume Scattering Function of 
ocean waters in the backward direction is dependent upon the wavelength (X). 

20 In a further embodiment of the present invention, the Volume Scattering Function of 
ocean waters in the backward direction is obtained at a particular wavelength using 
appropriate wavelength depenedent products of remote sensed reflectance and diffuse 
attenuation coefficient. 

25 In an embodment of the present invention, the Volume Scattering Function of ocean 
waters in the backward direction is obtained at the wavelength of about 490 nm. 

In another embodiment of the present invention, the method of interpreting the 
backscattered flux in terms of the VSF at a fixed wavelength is extendable to all the eight 
30 spectral bands of the Oceansat sensor which are located at 410nm, 443 nm, 490 nm, 510 
nm, 555 nm, 670 nm, 785 nm and 865 nm. 

In yet another embodiment of the present invention, the method can be extended to 
contiguous scan lines along the satellite track to generate a high resolution two 
35 dimensional volume scattering surface at fixed backscattering angles. 
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In still another embodiment of the present invention, the method can be extended to 
create a set of two dimensional volume scattering surfaces in several different visible 
bands by the satellite ocean color sensor. 

5 Accordingly the present invention provides a method of determining the Volume 
Scattering Function of Ocean Waters in the backward direction using a satellite Ocean 
Color Sensor. The said method exploits the geometry of the sun-ocean-satellite detector 
to function as a backscatter system comprising of the sun [1] as a source of radiation, the 
ocean as the sampling volume [2], and the satellite ocean color sensor [3] as the detector 

10 of the backscattered flux from the ocean emanating from separate ocean pixels arranged 
on a scan line across the track of the satellite payload and arriving at corresponding 
pixels on the CCD array detector of the satellite sensor [3]. The Applicants have 
surprisingly found that the output signals from each electronic pixel of the Ocean Color 
Sensor is related to the Volume Scattering Function at a fixed scattering angles which 

15 bears a direct relation to the cosines of the solar incident angle (G s ) and the angle (0 P ) of 
the upwelled flux at the ocean surface, the Remotely Sensed Reflectance (Rrs), and the 
diffuse attenuation coefficient (Kd) thereby resulting in a new product of the Volume 
Scattering Function at fixed backscattering angles. 

20 The method when extended, can result in a more complete spectral characterization of the 
Volume Scattering Function at eight fixed wavelengths and at 3944 fixed scattering 
angles determined by the view angles of the detector pixels and the incident solar angle. 

The invention is further described in the following examples. The Exampleas are 
25 provided to illustrate the present invention and hence, the same should not be construed 
to limit the scope of the present invention in any manner, 

EXAMPLE 1 : CALCULATING THE SHARPNESS OF THE FIELD OF VIEW 
(FOV) CONE ANGLE SUBTENDED BY AN OCEAN PIXEL AT THE 
30 SATELLITE DETECTOR PIXELS 

This one-time calculation demonstrates that the FOV subtended by a pixel on the ocean 

surface is extremely narrow, thus simplifying the method of computing the backward 

VSF by space sensors. Using the parameters in Table 1, we assume an ocean pixel of 

35 area 360m x 236m and a satellite altitude of 720 kms. The solid angle subtended by an 
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off-centre ocean pixel of area 360m x 236 m at the satellite detector height H ( = 720 
kms) is given by : 

dQ - [ 360x236 cos0 p ]/[H 2 ] ....[1] 

5 Calculation shows that the solid angles are of the order of <= 2.5 mrads for look angles 
0 P in the range +43 ° to -43 0 . These are extremely narrow FOV angles, and it is safe to 
assume that the pixel cone can be replaced by a straight line from the centroid of the 
ocean pixel to the CCD pixel on the satellite detector. 

10 EXAMPLE 2: CALCULATING THE SCATTERING ANGLE if/ OF THE 

BACKWARD VOLUME SCATTERING FUNCTION p(v|/, X) AT A GIVEN 

WAVELENGTH. 

The OCM sensor on the IRS-P4 satellite consists of eight CCD line array detectors, each 
15 of which view the backscattered radiation in the eight wavebands shown in Table 1. At 
any given instant of time, the CCD pixels on a detector view a scan line of ocean pixels 
arranged along a 1420 km swath of the ocean surface as shown in Fig. 3. Each of the 
3730 CCD detector pixels is associated with a look angle 6 P , which have values lying in 
the range + 43 ° of the fan shaped cone shown in Fig.3. Simple geometric anaysis of the 
20 representation in Fig.4 shows that the scattering angle \j/ can be defined as the angle 
between the extended solar ray direction and the direction of the backscattered water 
leaving radiance emerging at the look angle 0 P to the zenith direction by : 
V = [ 18O-(0 P +e s ) ] ....[2] 

25 Given the simple geometric relation in [2], and assuming 0 S ~ 17° as a typical solar angle 
shows that \\f will vary over the range 127° to 154° for look angles 0 P = + 43° of the 
total detector FOV. The solar angle 0 S changes with time of day and with Julian day, and 
therefore is a precise and predictable natural variable that controls the range of scattering 
angles \\f through the relation [2]. 

30 

EXAMPLE 3: CORRECTING THE LOOK ANGLE 0 P OF A CCD PIXEL ON AN 

OCEAN PIXEL 

The look angle 0 P of a CCD pixel on the satellite sensor is measured with respect to the 
35 zenith direction. However, a correction needs to be applied to this parameter by the 
relation ( see A. Morel, J.L. Mueller, " Normalised Water-Leaving Radiance and Remote 
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Sensing Reflectance : Bi-directional Reflectance and other factors " in Ocean Optics 
Protocols for satellite Ocean Color Sensor Validation, Rev. 3, Vol 2, pp 183-210, 2002) 
9 P = 0 S at +sin~ 1 (0.113 tan 0 sat ) ...[3] 

5 The angle 0 sat refers to the satellite view angle, and this is available in most satellite data 
processing packages for Ocean Color sensors 

EXAMPLE 4: PROVIDING AN OPERATIONAL FORMULA USED IN 
CALCULATING THE VOLUME SCATTERING FUNCTION P(\|/, X) IN THE 
10 BACKWARD DIRECTION AT A GIVEN WAVELENGTH. 

The derivation of the Volume Scattering Function P(\|/, k) is based on the definition of 
water leaving radiance L(0 P ) which forms the only component of backscattered radiant 
flux emanting from the ocean surface per unit solid angle per unit projected area. It can 
15 be shown that the VSF at a given wavelength is given by the product of the Remote 
Sensed Reflectance Rrs ( 0 P t A,), the diffuse attenuation coefficient Kd, and cosine proauct 
of the solar and look angles ( 0 S ,0 P ) respectively Thus 

P(\|/,X) = Rrs ( 0 p , A,) K<j [cos 0 s cos 0 p ] ....[4] 

20 The relation in [4] is independent of ocean pixel area, CCD pixel element or satellite 
altitude. Its units in m^.sr" 1 are that of a scattering function. Rrs is the ratio of the 
normalized water leaving radiance (nLw) to the mean extra-terrestial solar irradiance at 
mean earth-sun distance; Kd is the inverse of the satellite penetration depth (m); both 
quantities obtainable as derived satellite products. The spectral nature of Rrs and Kd 

25 makes P(i|/, X) a spectral variable also. The angular dependence of the VSF is controlled 
through the product [cos 0 s cos 0 p ] involving only the solar angle and the look angle. 
Given the the satellite view angle 0 sat , Rrs , Kd , and the solar angle 0 s at each ocean 
pixel on a scan line of the satellite swath, it is possible to use the VSF algorithm in [4], 
and the scattering angle in [2] to calculate a high resolution spatial distribution of the 

30 scattering function as a function of scattering angle along a scan line. For a complete 
swath width this provide 3730 points across the satellite track. The computation can be 
repeated for subsequent scan lines along the satellite track to reveal a complete 2D 
scattering distribution of the ocean surface. 
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EXAMPLE 5: APPLICATION OF THE VSF ALGORITHM IN [4] USING REAL 
SATELLITE DATA ACROSS A SCAN LINE AT 490NM WAVEBAND ( FIG.5) 

In order to illustrate the novelty of the algorithm in [4] , across track satellite pixel data 

5 was selected from an OCM image of 16 march 2002 when large patches of 

Trichodesmium blooms were detected in the coastal waters off the Goa Coast in the 

Eastern Arabian seaa. This particular dataset consists of all the quantitiesdetailed above 

and required in the computation of p(\|/, X) as afunction of the scattering vj/. The results 

are shown in Fig.5. Moving from the open ocean values of P(\|/, 490) ~ 0.0005 m^.sr" 1 at 

10 scattering angles = 138 degrees , the VSF shows a steadily increasing value until it 
encounters the Trichodesmium bloom at vj/ ~ 126 degrees where the scattering angle 
shows a significant increase to P(\j/ 5 490) > 0.0012 m^.sr" 1 near the coast. This simple 
example for a single scan can be extended for all contiguous scan lines along the satellite 
track , and for different wavebands so as to map the 2D spectral scattering distribution of 

15 the ocean surface. 

EXAMPLE 6: COMPARING IN-SITU BACKSCATTER DATA FROM A 
COMMERCIAL INSTRUMENT WITH THE ALGORITHM IN [4] USING 
SATELLITE DATA FROM THE SEAWIFS OCEAN COLOR SENSOR ( FIG.6) 

20 

The optical backscatter instrument described in the prior art above was used in a recent 
cruise to the Arabian Sea. This instrument measures the backscatter coefficient at a single 
backscatter angle. The upper plot in Fig. 6 shows the backscatter coefficient plotted for 
18 different stations of the cruise track. At station number 16 the ship encountered dense 
25 plumes of sediments in the water resulting in a large backscatter value of ~ 0.1 m" 1 . 

An independent approach based only satellite data from the SeaWifs Ocean Color sensor 
was applied using the algorithm of the present invention. Unlike the technology of the 
OCM sensor on IRS-P4, the SeaWiFs sensor has reduced ocean pixel resolution of 1km x 
30 1km, and does not employ the push broom technique used in OCM, At each station, 
SeaWiFs data on Rrs, Kd, 0 s 9 0 P was processed and the VSF algorithm applied to 
provide the lower plot in Fig.6 . At station 16, the VSF algorithm registered a ten fold 
increase in the backward volume scattering function exhibiting a significant sensitivity to 
the presence of sediments in turbid waters. 
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Advantages of the Present Invention: 

The main advantages of the present invention are : 

> Provides for a method of determining the backscattered portion of the Volume 
Scattering Function of oceanic waters by use of a satellite ocean color sensor. 

> Redefines the concept of a sample volume as that of an ocean pixel 360 meters wide 
and 250 meters long, and by doing so captures a true representation of suspended 
particles along large tracts of the ocean. 

> Extends the method of finding the VSF to eight other spectral wavelengths, and by 
doing so provides for additional information on the spectral nature of the VSF at 
ocean scales. 

> The computation can be repeated for subsequent scan, lines along the satellite track 
to reveal a complete 2D scattering distribution of the ocean surface. 
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